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Abstract 


The second order formalism for fermions provides a description of fermions that is 
very similar to that of scalars. We demonstrate that this second order formalism is 
equivalent to the standard Dirac formalism. We do so in terms of the conventional 
fermionic Feynman rules. The second order formalism has previously proven useful 
for the computation of fermion loops, here we describe how the corresponding rules 
can be applied to the calculation of amplitudes involving external fermions, including 
tree-level processes and processes with more than one fermion line. We comment on 
the supersymmetric identities relating fermions and scalars and the associated simpli¬ 
fications to perturbative calculations that are then more transparent. 
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The second order formalism for fermions provides a description of fermions that is very 
similar to that of scalars. The prefix second being due to the order of the derivative in the 
associated effective action—the ~ 1/ k 2 form of the propagator. Its has been obtained from 
an analysis of the infinite string tension limit of string theories BBS and is related to the 
first quantized form ||, [5|. It has also been used in perturbative calculation [j|], where it 
simplified the calculation of one loop amplitudes in QCD, gravity and the electroweak sector 
of the Standard model. 

It is the purpose of this letter to make clear the relationship between the first and second 
order fermion rules, in a manner that will be immediately obvious to Feynman diagram 
calculators. Further we generalize the fermion-loop inspired rules to include a description of 
external fermions. 

By way of an example, we shed some light on the supersymmetric similarity of scalars 
and fermions interacting with gauge fields. With respect to the mechanics of calculating 
Feynman diagrams this has previously been quite mysterious. 

For completeness, we begin by reviewing the formal approach to the derivation of the 
second order formalism. The basic idea being to exploit the relationship between fermionic 
actions and straightforward determinants. This is completely analogous to the derivation of 
the Faddeev-Popov ghost effective Lagrangian term 0, which is required for the calculation 
of general non-abelian amplitudes. 

We choose the following convention for our gauge coupling, 

D,j, = <9,, + igQA M . (1) 

For QED the coupling is g — e, and the associated charge is Q = —1 for the electron. In 
QCD, and related non-abelian theories, the coupling is required to be universal, Q = 1, 
and the gauge held is expanded in terms of hermitian generators for the gauge symmetry, 

= T a Af. With these conventions we define the held strength tensor to have the form, 

[D„ D v \ = (d,A„) - (d u AJ + igQ [A„, A v ]. (2) 


The starting point is the generating functional for the effective action [|[] of a closed 
fermion loop interacting with a gauge held, 


r [A\ = io gj 

1 ThfThf exp 

i J dx ip (ilf) — rnj if 

= log 

det [it> — m^j 

? 


( 3 ) 

( 4 ) 


where we have used the familiar Faddeev-Popov trick in reverse. Now, following the asso¬ 
ciative property of determinants and the anti-commuting properties of y 5 (det( 7 5 ) = 1), we 
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Figure 1: The basic fermion unit. 


can manipulate the above expression into the second order form 0 , 

1 


r [A] = ^ (log (det ip - m ) + log (det 75 (it* - ml) 75 )) 
- log det (ip — ml) (—ip — rnj 75 


- log det (p 2 + m‘ 


= ) log det ( D 2 + m 2 + ^-a^ u F, 


11V 


(5) 

( 6 ) 

(7) 

( 8 ) 


Here we use a^ v = |[ 7 ;i , A\- F° r convenience, we have simply reproduced the derivation of 
Strassler [0. 


Despite the factor of 1/2, this generating function can be thought of as describing a field 
theory for fermions. Not the original fermions, "0 of Eq. (||). Instead, a new set of fermions, 
T, whose effective action is, 


r [A\ = - log J DTDT 1 ’ exp 


1 1 dx vlT [d 2 + 


m 2 + —a^F, 
2 


/LV 




(9) 


The Feynman rules for such a theory can be derived from this exponent. Following an 
alternative derivation, we give them diagrammatically below, see Fig. [|. 


In order that we need not consider the external properties of the new fermions, T, we 
shall concentrate here on calculations of amplitudes involving an internal fermion loop. This 
closed loop may in general be embedded in a multi-loop diagram. We shall deal with external 
fermions afterwards. By reorganizing the conventional calculation for such an amplitude, we 
shall show how it can be equivalently described in terms of rules of the second order form. 
This approach will enable us to apply the second order rules to cases where the determinant 
formula does not apply. 


It is the case that all fermion-gauge interactions lead to the need for evaluating the trace 
over some spin line. Where the spin line is built from products of elements of the form given 
in Fig. [Tj. This is simply the product of a fermion propagator and a fermion-gauge vertex. 


2 













In an abelian theory such as QED, the factor T a is the charge of the fermion with respect to 
the gauge held, Q. For non-abclian theories, such as QCD, it corresponds to the appropriate 
generator for the coupled gauge held. The mass of the fermion is m and the gauge boson 
carries the momentum, —p, away from the vertex. 


It is our intention to manipulate this unit in such a way that we arrive at a set of 
equivalent Feynman rules whose principal property is that all of the Lorentz structure is 
contained in the vertices. For brevity we write this term as A k ’^/D k + P , where, 

Aj*’“ = igT a (f +f + m) 7 " and D k+P = i((k + p ) 2 - m 2 ). (10) 

By applying a Gordon-type manipulation j3j we can write A/'/ in the following form, 


At; = 'j-T‘{(2t + 2m + f)r+p a ) 

= igT° ((k“ + p”) - (-*:“) - 1[7“,)*] - 7° (t ~ ™)) . 


( 11 ) 

( 12 ) 


and accordingly we dehne, 

= igT a ((k a +p a )-(-k a ) + icT af3 py) and C k ’ a = igT a ^ a (ft - rnj , (13) 

such that, 


a a,a _ 


fycx.,ci 

L, k ■ 


(14) 


We observe that the form of C is exactly sufficient to ensure that it generates a propagator 
when followed by an appropriate A: 


£iQ!.,CL 


A^l q _ q = -1 x i(k 2 - m 2 ) x -ig 2 T a T b 7 “ 7 /3 = -D k x E a<i - 


ab 


(15) 


where, 

E *P,ab = _ ig 2 ^af3 _ m a.3\ 

These expressions serve to dehne E a ^ ,ab , which (after having removed the two color matrices, 
T a T b ) will correspond to a color ordered four-point fermion-gauge vertex. This color ordering 
arises from the fact that individual Dirac fermion loop diagrams contribute ordered color 
traces to the amplitude. It has been established elsewhere 0, that color ordering is a 
remarkably efficient tool for use in the calculation of non-abelian gauge theory amplitudes. 
To include abelian theories (QED), however, we shall direct our attention to the sum over 
permuted diagrams. 


We consider some general fermionic loops involving strings of these A/D objects. By 
re-writing selective As as (B — C ) we will generate an equivalent description of the fermion 
calculation in terms of a new set of Feynman rules. These rules are identical to those 
anticipated from the second order form of Eq. (|). 
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Figure 2: A generic fermion loop, as described in Eq. (|T7|). We label the legs from 0, 
such that the loop momentum flowing into the 0 -vertex along the fermion line is k, and 
Po = -£?=i Pi- 


The generic {n + l)-external boson fermion loop, as derived with conventional Feynman 
rules, can be written as a sum over the n\ orderings of boson couplings to the loop. The first 
of these orderings has the following form, 


A /j 

rr _ rp„ ) k+...pi,p n ^ l fc+...pi,p rl _2 

±n — 11 ' n n n 

. ±y k+pn...+pi ±y k+p n -l-..-¥pi - Ly k+p n -2...+pi 


fciPn-"iPl 5 


D L 


(17) 


where we have neglected the integration over the loop momentum. The minus sign is the 
familiar fermion loop factor. Replacing the first of the As in this trace with Eq. (0) leads 
to, 


T = 

- 1 n 


-Tr 


JD&n i&n 

n k + ...p\,pn 


D k 


Aa„-l,a n -i pu„o n -i,a„a„-i ' 

^k+.- .piVn-l & 


A (-%n — 2 j tin — 2 
/i fc+...pi,p rl 2 


/t«o,ao 

^ 1 fc+p n ...+pi,-p n ...-pi 


fe+p n ...+Pl -^fc+Pn-l...+Pl 


D 


k+pn...+p 1 


D 


k+p n -2---+Pl 


D, 


(18) 


This pattern of replacement can be repeated: working our way from left to right we substitute 
for the left-most A with (B — C) and then apply Eq. (|15|). What remains once we reach 
the final A is a series of products of B s and Es plus a different series of B s and E s that are 
truncated with a ...+ P1 - We note that all of the latter terms carry a leading minus sign. 
For example, in the n — 1 (self-energy) case we obtain, 

B?T KTpl-p k + E^^o _ BIT 

Dk+pi Dk Dk+p! Dk+p! 
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It is desirable to get rid of terms that contain the object C. And to do this we appeal to 
the cyclic property of the trace. We shift the trailing C of the offending term to the left and 
then apply the rule that any B immediately to the right of a C can be replaced by (A + C). 
Here this corresponds to, 


Tr 


>^<ao,ao jDOii,ai 
^ k-\-pi *^k,p\ 


Di, D 


= -Tr 


k+pi 


^jaiao,aiao 

Du 


(~<oto,ao fyai,a i 


Du D 


k+pi 


( 20 ) 


where again we have applied Eq. (|15D . By comparing Eqs. (ED and (0), it is simple to see 
that the cyclic property of the trace guarantees that, 


Tr {c Q C n C n - l C n ~ 2 ...} = Tr [A n A n ~ l A n - 2 ... A 0 }, 
which means that the —CC term, of Eq. (^), is exactly — Tj and thus we have that, 


Tj = —Tr 
2 


j^OLi,ai jgao,ao jjjcxicxo ,aiao ^QiQo,aiao 

k,p i k-\-p i, p i ^ ^ 


D 


k+p i 


Du 


D 


k+p i 


Lh 


( 21 ) 


( 22 ) 


The factor of 1/2 is quite general—for arbitrary n. For this simple example the last two 
terms are actually identical (this follows from the invariance of the appropriate loop integral 
under the change of loop momenta, k —> — k — p±, and the cyclic property of the trace). 
Quite generally, however, we can write the sum of two E s with their indices reversed as, 


Qa/3,ab __ j^ot\a 0 ,a\a 0 _|_ j^aia 0) aia 0 __ _^2 j r J' a _ 


%<J 


a(3 


rjna rjnl 


)■ 


(23) 


Now, since we are attempting to create a new set of Feynman rules we are at liberty to 
mix contributions from different first order Feynman diagrams. To this end we consider the 
joint contribution of Figs. 0« and |2jb. Amongst the many B and E terms that arise from 
the replacement of As, there are a series of terms that are identical except for the ordering 
of the indices on a single E object. These terms correspond to the relative reversal of two 
coupled bosons, as is depicted in the figures. 

It follows that any physical process described by a trace over terms of the form A/D , 
can equally well be described by a trace over terms built from the objects, T>, D and G. 
These new rules for Feynman diagram calculation are given in Fig. ||. We note that in the 
simplifying case of QED, T a —> Q, the four-point vertex reduces to, 

Gqed = - 2 %«)V' 3 . ( 24 ) 

We may immediately identify these rules with those derivable from the second order deter¬ 
minant, Eq. (HI). We note that the direction of the T fermion arrow is identical to that of 
the first order fermion, ip. That is to say, we label the head of the fermion line by Th It is 
also worth noting that the mass of the fermion is now only present in the denominator of 
the propagator and it no longer complicates the Lorentz structure of the Feynman rules. 

The reorganization of Feynman rules described above is crucially dependent on the fact 
that we can write the first order rules for individual diagrams as a product of the unit-terms 
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Figure 3: The second order rules for calculations of fermion-gauge interactions. Care should 
be taken when using these rules to ensure that the trace is taken over all terms and not just 
those containing a a a/3 . All momenta are directed outwards, which accounts for the sign 
difference between the three point a a/3 contribution and that of the B in Eq. (|T3|). When 
one traces over a closed fermion loop one should also multiply by 1/2. For internal closed 
fermion loops there is the familiar additional factor of —1. See the text for the application 
of these rules to processes involving external fermions. We note that by setting a a/3 —> 0 and 
i —> —i, these rules become those for a gauged-complex scalar field. 
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in Fig. [T] (As and D s). It is not immediately apparent that this is the case for processes 
with external fermions, where the associated amplitudes are constructed from one fewer 
propagators than vertices, sandwiched between a pair of spinors. That is to say, the leading 
spinor is immediately followed by a first order vertex and not a propagator, so the string of 
objects sandwiched between the spinors may not immediately be rewritten as a string of As 
and D s. 

At the level of the amplitude squared we should have confidence that the process may 
be described with respect to the second order rules: the loop contribution and the squared 
amplitude are directly related by the optical theorem. Indeed at this level, we create a trace 
by moving the trailing u (or v ) spinor of the amplitude to precede the leading u (or v ) 
spinor of the conjugate amplitude. Summing over the spins of the fermions we replace the 
associated spinor products as follows, 

'^2ui(p,m)ui(p,m) = f + m and ^ w;(—p, m)bj(—p, m) = — + m) (25) 

i i 

(where we give the argument momenta of the spinors relative to the direction of the arrow 
on the fermion line). These energy projection matrices are formed in exactly the correct 
positions to construct A objects from the problematic first order vertices; at the head of the 
amplitude and conjugate amplitude. 

The absence of a propagator in the denominator, 1/D, at the point in the trace where 
the external fermions are represented ensures that no Es are constructed across the cut. 
Hence, this reallocation of spinors is conveniently factorizable into two parts that we may 
choose to label as a second order amplitude and conjugate amplitude. Unfortunately, this 
procedure leads to a representation for the amplitude that is a 4 x 4 matrix, as compared 
to the more efficient complex number of the conventional approach. Further we have lost 
all of the spinor information for the fermions. A more useful approach, which we shall now 
describe, is to effectively insert the appropriate energy projection matrix after the leading 
spinor and so retain a simple scalar representation for the amplitude. 

By inserting an energy projection matrix, we can write any amplitude involving massive 
first order fermions in the following way, 

uAk, m )T ... = 1 (j/ + m) T ... (26) 

2m v 7 

(Note, antifermion spinors, v 3 (p, rn ), have momentum arguments that flow against the di¬ 
rection of the fermion arrow so in terms of the As they behave identically to the fermions 
but with an overall factor of —1 attached.) It follows, in terms of the objects introduced 
above, that each amplitude takes the form, — uA(A/D).. ,u/2m. Equivalently, this may be 
written as a spinor sandwich where the filling is constructed from the second order rules 
of Fig. The factor of 1/2 encountered when we dealt with closed loops is not required 
here because the trailing Uj(k,m ) annihilates any C*, that precedes it—cf. the discussion 
around Eq. (|T|). Operationally, we generate the factor of 1/2 by only using the fermion or 
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antifermion contribution to the field at the ends of the fermion line (depending on whether 
it is in the initial or final state). 


Much simplification in the calculation of held theory amplitudes has been achieved with 
the use of spinor-helicity methods [[ 


m m 


ft is desirable to see if and how such 
techniques may be applied with respect to the second order formalism for massless fermions. 


The trivial rewriting above for the massive fermion case collapses in the massless limit, 



following form, 



and once again we may re-express such an amplitude in terms of the rules in Fig. ^sandwiched 
between the arbitrarily chosen x-spinor and the original trailing p^,-spinor. The choice of x M is 
as arbitrary as that made for each of the reference momenta used to construct polarization 
vectors. In other words, it is only required to be non-coincident to the original spinor’s 
momentum, pA Note, as with the polarization vectors, we cannot choose different xs for 
each of the separate Feynman diagrams since we have mixed diagrammatic contributions to 
obtain the second order rules. 


An obvious choice for x^ is that of the trailing spinor, here given as p (0 As written 
out in Fig. [5[ the vertex rules are in general the sum of two distinct pieces: a piece which 
is essentially identical to the rule for a simple scalar interaction and a uniquely fermionic 
Lorentz -structure carrying term, ~ cd"', see Fig. |j. The choice x = p^, in a single stroke, has 
the effect of cancelling the purely scalar contribution to the fermionic amplitude, because 
I <? A ) = 0. Below we shall make a different choice to account for a supersymmetric identity 
at a diagrammatic level. 


Since the second order rules bear such a close resemblance to the Feynman rules for a 
scalar interacting with a gauge held, it is worth investigating whether they shed light on the 
the supersymmetric Ward identities |T|, |T^]. We shall, for the sake of example, consider the 
simple relationship between the A( —f; + amplitudes of the 0^0 + (n — 2)7 and 

00 + (n — 2)7 processes at tree level. Here the 7 s represent abelian photons, 0 a complex 
scalar and 0 a fermion. We use the convention that the positive (negative) helicity label for 
the scalar refers to its particle (anti-particle) state. The non-abclian amplitudes obey the 
same identity [O], but for the sake of simplicity we shall only present the abelian case. 


We shall need the following identity, 

( X (p;r),f> 


4A 

V2 


-A 


P 


p 


(28) 
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Figure 4: This figure represents contributions to the MHV amplitudes which enjoy a close 
supersymmetric similarity. The text describes the simple relationship between the (a) second 
order fermionic and (b) scalar amplitudes. 


Here we use the polarization vector |T2], 113 


«i(p; r ) = A - 


(p x 

7m 

r A ^> 

V2( 

rp — A 

p x y 


(29) 


where p is the physical momentum of the photon carrying helicity A, and r is the reference 
momentum. 


Now, with reference to Fig. [|, we shall compare the product of terms used in each of 
the diagrams that contribute to the fermion (Fig. [|a) and the scalar (Fig. (4]b) amplitudes. 
Such amplitudes are termed maximally helicity violating (MHV), because they are those 
amplitudes with the most external legs of the same helicity that are not identically zero at 
tree level. 


In terms of the second order rules, we can see quite quickly that the M^(++;... +) and 
A^( —b;... +) amplitudes vanish. The former because each of the rules is helicity conserving 
and p+'j — 0 (this result is also obvious when using the first order rules). The latter 
follows from the following fact, 


n 


[{*{p\r),f} = 0 


(30) 


(see Eq. (^))— that is to say, where all the photon helicities are opposite to that of the leading 
spinor the fermion and scalar vertices are effectively identical. Choosing the arbitrary spinor 
momentum at the head of the fermion line to be equal to the momentum of the trailing 
spinor, x = gives an amplitude proportional to (p^\ p^j / (p^\ p^j = 0. In this latter 


case the proportionality factor is the scalar amplitude, A^( —f;...+). This can be seen 
to be zero by choosing the reference momenta of every polarization vector to be equal to 
momentum of the trailing spinor, pp. So in fact both of these fermion amplitudes are zero 
independent of our choice for x. This can also be seen as a consequence of supersymmetric 
Ward identities 
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For the MHV amplitudes the relationship between the scalar and fermionic amplitudes 
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is almost as simple as that of the A{ —f;... +) case. This is because each of the three-point 
vertices that couple the positive helicity photons to the fermion are effectively just scalar 
vertices (see Eq. ([30|)). The difference between the fermion and scalar amplitudes is due 
solely to the uniquely fermionic component of the negative helicity three-point vertex. This 
contribution is of the form (see the first of the graphs Fig. [3]a and Eq. (ffij|) ), 


x~ 


W / - 
1 „+\ \ x 

^before 

A>) A-> 

‘S'after 






Pi, 


(31) 


where p(_) is the momentum of the negative helicity photon. Sbefore (after) is the product of 
the preceding (trailing) scalar-like vertices and propagators, which freely commute with (x _ |. 
Clearly the choice, x = P(~), removes this contribution from the amplitude. 


Using Eq. (j27|) with the choice n = £>(_) for the spinor reference momentum, it becomes 
trivial with the second order rules to demonstrate in a perturbative sense that each of the 
diagrams contributing to these tree level MHV amplitudes (and accordingly the amplitudes 
themselves) are related by, 

•A*(-+;= j P( —>] (32) 

(*'<->! P*f) 

It is amusing to note that this diagram-by-diagram equivalence is even independent of the 
choice of reference momenta for the photon polarization vectors. At the level of the amplitude 
this is of course in complete agreement with the appropriate supersymmetric Ward identity 

ill. 


In this letter, we have discussed the second order formalism for fermions and have suc¬ 
ceeded in rewriting the Feynman rules for gauged-fermions in a form that exposes their close 
similarity to gauged-scalars. We have shown how these rules may be applied to the calcula¬ 
tion of both closed fermion loops and for processes involving external fermions. Because of 
the close relationship between the second order rules and the first quantized form in, this 
work may provide some reference for those wishing to construct string-like rules for external 
fermions. 


The supersymmetric relationship between fermion and scalar loops has been investigated 
previously i II where the second order formalism, as it applied to closed fermion loops, 
was employed. Here, we have demonstrated that this similarity exists for external fermions 
too. We have used it in a simple example, to make transparent the supersymmetric Ward 
identity relating two MHV amplitudes, which is something not normally obvious at the level 
of individual Feynman diagrams. 

Although the simple relationship presented in the above example is more complicated 
for other helicity configurations, it will always be the case that a fermion amplitude has a 
component which is directly proportional to a scalar one—a feature which is explicit with 
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the second order Feynman rules. This component is physical, insofar as it is gauge invariant. 
So by choosing a different set of reference momenta for the external gauge bosons, we may 
optimize the calculation of the fermionic remainder. That is to say, such a remainder may 
be calculated more efficiently on its own than when mixed up with the rest of the ampli¬ 
tude. This method of breaking-up amplitudes has already been used successfully for closed 
fermion loops 0, and can be thought of as a generalization of the established supersymmetric 
techniques 0.0. 
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